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ABSTRACT 

A 25 megabit/second direct detection optical communication system that 
used Q=4 PPM signaling was constructed and its performance measured under 
laboratory conditions. The system used a single mode AlGaAs laser diode 
transmitter and low noise silicon avalanche photodiode (APD) photodetector. 
Comparison of the measured performance of the system with that was theoreti- 
cally expected revealed that modeling the APD output as a Gaussian process 
under conditions of negligible background radiation and low (less than 10 -12 A) 
APD bulk leakage currents leads to substantial underestimates of optimal APD 
gain to use and overestimates of system bit error probability. A procedure is 
given to numerically compute system performancejwhich uses the more accurate 
Webb’s Approximation of the exact Conradi distribution for the APD output sig- 
nal that does not require excessive amounts of computer time (a few minutes of 
VAX 8600 CPU time per system operating point). Examples are given which 
illustrate the breakdown of the Gaussian approximation in assessing system per- 
formance. This system achieved a bit error probability of 10 -6 at a received sig- 
nal energy corresponding to an average of 60 absorbed photons/bit and optimal 
APD gain of 700. The measured performance of the system was found to be in 
excellent agreement with the performance predicted by the nearly exact computa- 
tional procedure used. 


* Work supported by the National Aeronautics and Space Administration 


I Introduction 

Recent advances in the technology of both semiconductor laser diodes [l]-[4] 
and low noise silicon avalanche photodetectors (APDs) [5]-[7] have generated con- 
siderable interest in free-space direct detection optical communication between 
satellites. Theoretical analyses of optical communication systems (both free-space 
and fiber-optic guide) that contain an APD as a photodetector have appeared 
previously [8]-[l0]. Due to the complicated nature of the statistical description of 
the response of the APD to incident photons, nearly all these analyses have 
modeled the APD output as a Gaussian process. The actual distribution of secon- 
dary electrons output by the APD in response to the absorption of primary pho- 
tons has been given by McIntyre and Conradi [ll] [12]. Use of the Gaussian 
approximation in the calculation of system bit error probabilities has been shown 
to be quite accurate in situations where the bulk leakage current of the APD 
itself is of order of nanoamperes, the APD gain is a few hundred, and background 
light levels are not negligible [13], [14]. In systems that use on-off keying as the 
modulation format, the optimal receiver consists of a threshold test [15]. In this 
instance, modeling the APD output as a mixture random variable consisting of 
the discrete number of APD output photoelectrons and continuous Gaussian ther- 
mal noise as an equivalent Gaussian random variable has been theoretically 
shown to be quite accurate, provided certain conditions are met [16]. 

The modulation format of choice in semiconductor laser diode direct detec- 
tion optical li nks consists of low order PPM signaling [17]. The optimal receiver 
[15] consists of a device which determines the time slot which contained the larg- 
est APD output signal; no comparisons against fixed thresholds are required. 
Furthermore, silicon APD’s with bulk leakage currents below 1.0 picoampere and 
values of k eff (ratio of the ionization coefficients of holes and electrons) of less 
than 0.01 are now readily available commercially [18]. These devices are conven- 
tional silicon APD structures and do not require fabrication of superlattice dev- 
ices in which k eff =0. 

A theoretical analysis of the measured performance characteristics of a 
AlGaAs laser diode direct detection optical communication system that contained 
a low noise APD photodetector will be described here. It shows that use of the 
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Gaussian approximation for the APD output can lead to substantial underesti- 
mates of optimal APD gain and substantial overestimates of receiver error proba- 
bilities under conditions of very low background light levels likely to be encoun- 
tered in satellite-to-satellite or deep space optical communication links. The rea- 
son for the discrepancy is that the Gaussian approximation does not accurately 
describe the low noise APD output signal in time slots that contain only low level 
background light. In addition, an efficient numerical computation procedure will 
be given for calculating receiver error probabilities that does not involve use of 
the Gaussian approximation. This procedure is shown to be in good agreement 
with measured system bit error rates. 

The remainder of this paper is organized as follows. The next section gives 
the details of the nearly exact numerical computation procedure for the determi- 
nation of bit error probabilities that does not require excessive amounts of com- 
puter time. It is shown that as the number of background light counts per time 
slot is decreased, use of the Gaussian approximation causes larger and larger 
discrepancies in predicted values of optimal APD gain and system bit error pro- 
babilities. The last section describes the prototype laser diode optical communi- 
cation system which was operated at a source data rate of 25 megabits/second, 
used Q=4 PPM signaling and achieved a bit error probability of 1CT 6 at a 
„ received signal energy corresponding to, on average, 60 absorbed photons per 
information bit (120 photons/PPM symbol). The theoretically computed perfor- 
mance is shown to be in good agreement with the experimentally measured per- 
formance. The nearly exact procedure for calculating system performance 
correctly predicts not only the optimal APD gain (which can be found experimen- 
tally), but actual system error probabilities as well. 


II Theoretical Analysis 

The probability that n photons are absorbed from an incident optical field of 

known intensity, P 0 (t) watts, by a photoelectric effect detector over an interval 

t+T 

[t,t+T] is given by a Poisson distribution of mean n=-|j / P 0 (t' )dt' , where 77 

“ t 

is the quantum efficiency with which the detector converts incident photons to 
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photoelectrons, and hf is the incident photon energy. In the case of APD photo- 
detectors, primary photoelectron-hole pairs undergo an avalanche multiplication 
process which results in the output of m electrons from the APD in response to 
the absorption of, on average, n primary photons. The integer random variable m 
cannot be less than the number of absorbed photons and is characterized by the 
Conradi distribution [12] 
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where G is the average APD gain. Due to the randomness of the gain mechanism 
in the APD, the output signal is characterized by an excess noise factor defined 
as F=E{m 2 }/E 2 {m}, which can be shown to be [12] 

F = k et G+(2-l/G)(l-k, a ) (2) 

Silicon APDs have quantum efficiencies of nearly 80% at wavelengths around 800 
nm, and values of k eff as low as 0.006 [7]. In addition, all APD’s have both bulk 
and surface leakage currents which are nonzero. This results in nonzero output 
current in the total absence of incident photons. Currently available low noise sil- 
icon APD’s have surface leakage currents of about 10 nanoamperes, and bulk 
leakage currents of 1.0 picoamperes or less [18]. 

In Q-ary PPM signaling, L binary source bits are transmitted as a single 
light pulse in one of Q= 2 L possible time slots once every T seconds. The output 
waveform has constant average power and peak power Q times larger. Semicon- 
ductor laser diodes are both peak and average power limited. As long as Q is not 
so large as to cause facet damage in the laser, PPM signaling produces high peak 
power short duration pulses that become more easily distinguished from back- 
ground and thermal noise as Q is increased. Peak power limitations of the sem- 
iconductor laser generally limit Q to values of 8 or less [17] [19]. 
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The optimal receiver for this communication format consists of a device 
which integrates the APD output signal over each time slot and then chooses the 
largest as the slot that contained the received light pulse. Figure 1 shows a block 
diagram of the system. The PPM word error probability (WEP) can be written 
as [9] 

00 x 

WEP = 1 - / p(x I X,T S ) [ / p(x' I X 0 T s )dx' ]«-■<& (3) 

-oo — oo 

and the corresponding bit error probability (BEP) as 

WEP (4) 

p(x | XT S ) is the probability density function for the integrated output from each 
time slot. represents the photon absorption rate due to received signal and 
background radiation and X 0 is the photon absorption rate due to background 
radiation alone. T s is the slot time period. 

The electrical charge integrated over each time slot is defined as x. Condi- 
tioned on the number of secondary photoelectrons, p(x | XT S ) can be written as 

00 

p(x | XT S ) = p(x | m)P(m | XT S ) (5) 

m=0 

The accumulated charge during each slot consists of the APD output electrons * 
and any additive amplifier thermal noise charge. Given the secondary electrons, 
m, x is a Gaussian random variable with probability density 

(x-X m ) 2 

p(x I m) = <£(x,x m ,cr 2 ) = - ^i-^ e 2<7 * (6a) 

x m = me + I S T S (6b) 

r = (2el* + )BT S 2 (6c) 

In equation (6), e = electron charge, K = Boltzmann’s constant, T = equivalent 
noise temperature, I s = APD surface leakage current, and B = one sided noise 
bandwidth. The APD surface leakage current, Ig, dose not get multiplied by the 
APD gain and therefore can be treated as a constant DC current. The APD bulk 
leakage current does get multiplied by the APD gain and can be modeled as part 
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of the background radiation. 

P(m | XT S ), given by (1) with n = \T S , has been approximated with good 
accuracy by Webb [20] as 


P(m | XT S ) = 
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for m values greater than XT S . Equation (7) approaches a Gaussian density with 
mean GXT S and variance FG 2 XT S for values of m close to its mean, i.e., 
| (m-GXT s ) | «GXT S . Figure 2 plots (7) and its corresponding Gaussian 
approximation, for XT S equal to unity and 120 respectively. As can be seen, the 
true distribution of secondary electrons skews toward the right and departs 
greatly from a Gaussian at both tails. These are the regions which form the 
main contribution to the PPM word error probability. It is also noticed that at 
very small mean number of primary incident photons, e.g. XT S =1, The Gaussian 
shape is not even close to the shape of the true distribution. Even worse, the 
Gaussian extends into the region of negative numbers of secondary electrons, 
which is clearly unphysical. Nevertheless, due to its simplicity, the Gaussian 
model is widely used in evaluating error probabilities of optical communication 
systems. 


If the distribution of secondary electrons, m, is assumed Gaussian, the 
charge accumulated during each time slot, which is the sum of secondary elec- 
trons and amplifier thermal noise, is also a Gaussian random variable. Therefore, 
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where X is given by \+X b +I b /e or X b +I b /e. X b represents the photon absorption 
rate due to the actual background light and I b /e represents the contribution of 
the APD bulk leakage current to the APD output. Here B is taken to be l/2T s 
and is the equivalent noise bandwidth of an ideal integrator. The power spectral 
density of the shot noise from the APD is taken to be d<i 2 >/df = 2e 2 XG 2 F as 
given by Smith and Personick [8]. 


It was found during the course of our experiments that at low background 
radiation levels, for example, (X b +I b / e)T s ~0.01, the use of the Gaussian approxi- 
mation significantly underestimated the optimal average APD gain and overes- 
timated the error probabilities. That can be understood as follows. In the 
extreme situation when (X b +I b /e)T s =0, the APD output from time slots where 
no light pulse was transmitted consists of only the amplifier thermal noise which 
is independent of APD gain. The error probability can then be reduced by 
increasing the APD gain, because the distribution of the output from those slots 
that do contain light pulses will be shifted and skewed to the right, away from 
the distribution of the output from slots that contain no light pulse. That 
implies an optimal APD gain of infinity and an optimal error probability which is 
infinitesimal. On the other hand, when the Gaussian approximation is used, 
increasing the APD gain not only shifts the distribution but also spreads it out 
symmetrically rather than skews- it towards the right. The probability that the 
APD output is smaller than a fixed threshold x 0 is 
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where 

T = x 0 -eGX s T s 

(e 2 G 2 FX,T s )‘/2 

Here the amplifier noise has been ignored. It can be seen that T— »-0 as G— >-oo for 
any value of x 0 because F is proportional to G, and consequently 
P(x<x 0 ) — **1/2 as G— *-oo. That is to say that the optimal APD gain cannot be 
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infinite and the minimal error probability never reduces to the infinitesimal. In 
reality, X b >0 and I b >0 but as long as (X b +I b /e)T 3 is small, the Gaussian approx- 
imation still underestimates the optimal gain and overestimates the error proba- 
bilities. The smaller (X b +I b /e)T s gets, the bigger the discrepancies become. 

Very low noise APDs have bulk leakage currents of less than 1CT 12 amperes 
[18], which is . equivalent to a background radiation count rate of order 
10 6 /second. Background light levels in space are negligible unless the sun or the 
earth enters the field of view of the receiver optics [21]. At useful data transmis- 
sion rates of 100 megabits/second and higher, the mean number of background 
counts per slot, (X b +I b /e)T s , can be well below unity. The Gaussian approxima- 
tion does not accurately model the APD output under these conditions. Conse- 
quently, a more accurate approximation to the Conradi distribution for the APD 
output must be used to correctly predict the system performance for very low 
background noise levels. 

Substituting (5) and (6) into (3), the PPM word error probability can be 
written as 
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x m and x m / are the mean values of charge x for given secondary electrons m and 


m' respectively. Since Erf( 


X-X m ' 


■) = 1-Erf( - 


x-x m / 


•), X m / «x m , and x 


a a 

varies about its mean during the course of the outer integration in (10), 
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)«1. We can then use a Talor expansion for the Q-l power and 


approximate the expression inside the braces of (10) and rewrite WEP as 
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00 


00 
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m=0 


-oo 


00 x-Xm/ 

• £ P(m ; I X 0 T s )Erf( — )dx (12) 

m' — 0 a 

This is equivalent to the union bound for the word error probability. 

In a numerical evaluation of (12), the two infinite sums and the infinite 
integral have to be truncated to finite terms and limits. Equation (12) can be 
rewritten as 
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The resultant error consists of three parts corresponding to each truncation. A 
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set of bounds for those errors are given in Appendix, along with the procedure for 
choosing the appropriate limits of the sums and the integral. The procedure to 
evaluate equation (13) is as follows. First, the limits of the integral for a fixed 
value of the index of the outer sum, m, Anj and B m , are computed. Corresponding 
to each value of the integration variable, x, the inner sum and then the integrand 
are evaluated. The integration is multiplied by P(m | X^) and the product is 
then accumulated. The process repeats for all values of m from zero to M. The 
subroutines for the integration and error function Erf(u) in (13) were called from 
the IMSL library. 

The use of Webb’s approximation, (7), for P(m j \T S ) in (13) rather than the 
exact Conradi distribution given by (l) greatly reduces the computational time 
needed to evaluate the WEP. However, Webb’s approximation is only valid for 
m>XT s >0 and terms corresponding to mf^Tg and m' <X 0 T S have to be 
treated separately. The contributions from the terms for m^XjTg in the outer 
sum of (13) are negligible since P(m | XjTg)— *0 for values of m which are much 
less than its mean, GXjTg. On the other hand, the contributions from the terms 
for m' =0 in the inner sum of (13) cannot be neglected when X 0 T S «1. The 
probability that no secondary electrons are output by the APD, P(m' =0 | X 0 T S ), 
is given by e 0 s , the same as the probability that no primary electron is gen- 
erated based on the Poisson distribution. If both X b T s and IbT s /e approach *2ero, 
which occurs at high data rates , low background radiation level, and low APD 
bulk leakage current, the probability that no secondary electron is output by the 
APD approaches unity. In other words, when X 0 T S «1, the contribution of the 
term for m' =0 in (13) must be taken into account and equation (13) should be 
modified as 
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where IgTg is the value of x m / at m' =0. P(m | X^) and P(m' | X 0 T S ) are given 
by Webb’s approximation, (7). The contribution from the second term in (14) 
increases as X 0 T S decreases. In the extreme case when X 0 T s =0, the second term 
in (14) becomes the only contribution to WEP and the first term vanishes 
because 

M M' X— X / 

E £ P(m' | X 0 T,)Erf( P-)dx 

m=X 1 T, A m m' =1 a 

oo 

< £ P(m' | X 0 T S ) 

m y =1 

= 1 - e " x ° Ts — <-0, as X 0 T s -+0 

Even more computation time can be saved by increasing the summation 
indices in (14) by increments bigger than unity with subsequent multiplication of 
individual terms by the value of the increments. Special care must be taken for 
small values of m' where P(m' | X 0 T S ) changes rapidly as m' increases. 
Experience showed that m' should be incremented by unity over the range 
1 <m r <100. Trial and error methods indicated that differences in the results 
from trial to trial could be maintained below 1 % for increments of m and m' 
(when m' >100) which were of the order of several hundrgd. 

Ill Experiments and Measurements 

A prototype optical communication system was constructed with the use of 
of a GaAlAs laser diode (Sharp LT024MD) as the transmitter and a low noise 
APD (RCA C30902S) as the photodetector. The system used Q=4 PPM signaling 
and operated with a source data rate of 25 megabits per second. A detailed 
description of the system is given in [22]. The integrations shown in Figure 1 
were realized with the use of a matched filter which had a tapped delay line 
structure described in [22]. Both the PPM modulator and demodulator shared 
the same timing signal (i.e. the same slot clock) for the sake of simplicity. A tim- 
ing recovery system must be provided at the receiver in an actual system. 
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The bit error probabilities of the system were measured as follows. A pseu- 
dorandom binary sequence 1023 bits long was input to the PPM modulator. The 
recovered binary bit stream from the PPM demodulator was then directly com- 
pared with the properly delayed version of the transmitted binary sequence with 
the use of an exclusive OR gate. The gate output was sampled at the source bit 
rate to avoid miscounting burst errors. A counter then read out the average error 
frequency. The bit error probability was computed by dividing the average error 
frequency by the bit rate. The counter recored at least 100 error events before an 
average error frequency was computed. 


The average optical power received by the APD was measured by removing 
the APD and substituting the sensor of an optical power meter which had a 
much larger active area than that of the APD. No other changes were made to 
the optics during the measurements. The APD was then repositioned and the 
average error frequency was again determined. The data were abandoned if the 
two average error frequencies had changed by more than 10%. This was done to 
avoid the effects of any drift in optical power during the course of the measure- 
ments. The average number of signal photons contained in each PPM word, 
which is the same as the average number of photons contained in each laser pulse 
if X 0 «\ 1 , was determined through the relation 


XiT s 


*?QPavT s 

hf 


(15) 


where P av is the average incident optical power in watts registered by the optical 
power meter and the quantum efficiency, 77 , was taken to be 77%. There are 
L=log 2 Q source bits in each PPM word and consequently the average number of 
photons for each bit is X 1 T s /log 2 Q. 

The optimal average APD gain for a fixed input optical power level was 
determined by adjusting the APD reverse bias voltage until the error probability 
reached a minimum. The actual value of the APD gain was measured as follows. 
The APD preamplifier output signal was displayed on an oscilloscope, and the 
average value of the peak height of the output in time slots that contained the 
received light pulse was estimated by inspection. Since the conversion factor, g 
(transimpedance amplifier gain in volts/amperes), of the APD preamplifier was 
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known, the average APD gain is given by 


G - — V /i 

eT ?QPav/hf 


( 16 ) 


where V is the average peak amplitude of the APD response in volts. 

Amplifier noise can be neglected compared with the APD excess noise if a 
strong cw optical field illuminates the APD. The excess noise factor, F, can then 
be determined by measuring the noise power output by the preamplifier over a 
known bandwidth. The noise power spectral density is given by [13]. 

= 2e -_G2F (17) 


Here P cw is the incident optical power measured with an optical power meter and 
d<i 2 >/df is determined by dividing the noise power by the bandwidth over 
which the noise power was measured. This method assumed the noise output 
was white, which was confirmed by direct measurement of the power spectrum of 
the APD output signal over the range dc to 100 MHz under conditions of con- 
stant intensity illumination. The excess noise factor F was computed from (17), 
and (2) was then used to compute k eff . Measurements made at several values of 
APD gain gave values of k eff of 0.010. 

The APD load resistor was 1030 fl as given by the manufacturer of the 
preamplifier and its effective noise temperature was found to be approximately 
1100 '’K. This result was obtained by measuring the noise power output of the 
preamplifier when the APD bias voltage was reduced to nearly zero (i.e. G^0) so 
that only amplifier thermal noise was present. 

Equation (17) was also used to determine the background radiation levels 
through the relation 



1 d_ , 2 

2e 2 G 2 F df >b 


(18) 


Here <i 2 > b represents the part of the noise power density contributed by the 
actual background radiation and the bulk leakage current of the APD. It is found 
as the difference between the total measured noise power density and the contri- 
bution from the preamplifier determined previously. However, (18) could not be 
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used to determine accurately the equivalent background noise count rate due to 
the bulk leakage current of the low noise APD, I b /e, by eliminating all the real 
background light (X b =0), because the noise at the APD output due to the APD 
bulk leakage current was only a few percent of the amplifier noise. However, I b /e 
could be calculated from the data supplied by the manufacturer of the APD for 
the spectral noise density in dark, i.e. (d < i 2 > / df) 1 / 2 = 3 X lCT^A/Hz 1 / 2 at 
G=600 [18]. I b /e is then given by the right hand side of (18) which was 0.0122 
per PPM slot time (20ns) corresponding to I b =9.79 X 10 _14 A. The APD surface 
leakage current is found by using the relationship 

Id = GI b + I s (19) 

where I d is the total APD dark current. For the APD used, I d was 12X10 -9 A as 
stated by the manufacturer, the APD surface leakage current was then 
11.9 X10 -9 A, about the same as the total dark current. The noise density due to 
this surface leakage current at the front end of the preamplifier was 
(2eI s ) 1 / 2 =6.2 X lO'^A/Hz 1 / 2 , and the noise density due to the APD load resistor 

i 

of 1030 O at the equivalent noise temperature 1100 °K was 
(4KT/R) 1 / 2 =7.3X10 -12 , much larger than that due to the surface leakage 
current. Therefore the APD surface leakage current had little effect on the meas- 
urements and system performance. 

Figure 3 shows the results of numerical evaluations of the bit error probabil- 
ity as a function of the average APD gain based on equation (14) for the optical 
inputs indicated in the graph. The curve with X 0 T S =12.5 corresponds to 
I b =0.1nA if X b =0, and the curve with X 0 T s =O.O122 corresponds to our actual 
experimental conditions where no real background light was present. The results 
based on the Gaussian approximation are also plotted (dashed curve) for the pur- 
pose of comparison. The measured optimal APD gain with X 0 T S =0.0122 was 
700, which was in good agreement with the numerically computed results based 
on (14) as shown in Figure 3. The Gaussian approximation works well for 
X S T S =12.5 but very poorly for X 0 T S =0.0122. The validity of equation (14) was 
further tested experimentally under the condition X 0 T S =1, which was set by arbi- 
trarily introducing background light according to (17). The measured optimal 
APD gain was 580, close to what the numerical evaluation had predicted, as 
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shown in Figure 4. 

Figures 5 through 7 are the numerical results of bit error probabilities as a 
function of the average number of photons required to transmit one bit of infor- 
mation, i.e. n s /bit= \ 1 T s /log 2 Q, under \ 0 T S =12.5, X 0 T 3 =1, and X o T s =0.0122 
respectively. Figure 5, once again, shows the appropriateness of Gaussian 
approximation under X 0 T S =12.5. The small crosses in Figure 6 and Figure 7 are 
the experimental results under each background noise level, and they are all very 
close to the numerical computation results based on equation (14). These results 
substantiate the validity of the numerical evaluation procedure described in the 
previous section. The Gaussian approximation is accurate only for high levels of 
background noise level and gets worse and worse as the average number of back- 
ground noise counts decreases. At X o T s =0.O122 and X 1 T S =110, use of the Gaus- 
sian model underestimated the optimal APD gain by almost a factor of two and 
overestimated the bit error probability by more than two orders of magnitude 
over what was actually measured in the experiments. 

In the numerical computations, the subroutine for the integration in (14) was 
set to have a relative error less than 1 %. The limits of the two sums in (14), M 
and M' , were chosen according to the procedure given in Appendix. As an exam- 
ple, at X 1 T S =110, X 0 T S =0.0122 and G=700, M~ 119,000 and M # ^133,000, the 
proper increments for the indices m and m' (when m 7 >100) of the sums in (14) 
were 2,200 and 500 respectively. Therefore about 50 points were taken to esti- 
mate the outer sum and about 350 points (including 100 points for 0<m 7 <99) 
were taken to estimate the inner sum. As a result, the CPU time required to 
compute one value of WEP was about two minutes on a VAX 8600 computer. 
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Appendix 

Limits of the Summations and Integration in (13) 

From Eq.(13), the truncation error is 

e = e x + e 2 + e 3 (Al.a) 

e i = (Q-l) E P( m I XiT s ) (Al.b) 

m=M+l 

oo — 

• f g P(m' I X 0 T s )Erf(-^l)dx 

-oo m' =0 a 

M 

% ■ (Q-l) E P(m | X,TJ (Al.c) 

m=0 

oo X— X / 

/ <A(x,Xj ^a 2 ) E P(m' | X 0 T s )Erf( — )dx 

x>B m> x<A m m' =0 ° 

M Bm 

e 3 = (Q-l) E p (“ I X,TJ / (Al.d) 

m=0 A m 

• E p (”' | X 0 T 5 )Erf(-^=E)dx 

m' =M' +1 ® 

where M and M 7 are the upper limits of the Sums, and A m and B m are the lower 
and upper limits of the integral. Since probabilities do not exceed unity and 
x-x m / 

Erf( )<1 for any value of x, 

<r 

M 

e 2 < (Q-l) g P(m | X^)- / <^(x,x m ,c7 2 )dx (A2.a) 

m=0 x>B m , x<A m 

The integrand in (A2.a) is now a Gaussian density function and A m and B m can 
be chosen such that the integration is a constant for any value of m. Then, 

e 2 < (Q-l) / <j>(x,x m ,cr 2 )dx (A2.b) 

x>B m , x<A m 

A m and B m can be set to be symmetric about x m and deviate by an amount 
which is determined by the specified maximum value for e 2 . For example, 
A m =x m -6cr and B m =x m +6a for e 2 <(Q-l)XlCT 9 . 
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Next, since x and x m i become further apart as m increases, the value of the 
integration in (Al.b) decreases monotonically as m increases. Therefore, 

oo 

e i < (Q-1) E P(“ I x iT,) 

m=M+l 

oo — 

oo X— X I 

• / ^(x^M+i^) E PK I x o T s)'Erf( — — )dx 

-oo m' =0 ° 

00 — 
oo X— X / 

<(Q— l)/0(x,x M+ i,a 2 ) g P(m' IVT S )-Erf( — )dx 

-oo m' =0 a 

Bmw M' X-X / 

«(Q-1) J ^(x,? M+1 , cr=) £ PK I X 0 T S )-Erf( l=-)dx (A3) 

Am+i m' =0 ° 

Hence it is safe to choose M such that the integral given in (A3) is less than a 
specified value (e.g. 10 -9 ). The integral in (A3) is identical to that in (13) when 
taking m=M+l, so proper values of M can be determined during the course of 
the numerical evaluation of (13), rather than through trials of separate computa- 
tions of (A3). Every time the index of the outer sum is increased, i.e. a new term 
is computed and accumulated for the outer sum of (13), the values of the integral 
within this new term is compared with a preset limit, for example, 10 -9 /(Q-l) for 
ejClCT 9 . The process terminates when the value of the integral in the next new 
term becomes less than th^ preset limit. 

x-x m / 

Finding the proper value of M' is lengthy. First, since Erf( )<1 for 

C 7 

any value of x, 

M Bm n oo 

e3 < (Q-1) E p ( m I >'iT s )/#(x,x nl ,<T 2 )dx- Y PK IVJ 

m=0 A m m J =M! +1 

oo 

< (Q-1) Y PK I *oT s ). (A4.a) 

m / =M / +1 

Substituting (7) for P(m' | X 0 T S ) in (A4.a) and then bounding the sum by the 
associated integral gives, 
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00 


e 3 < (0-1) s 


M' (2rX 0 T s G 2 F)^ 2 m' -GX 0 T, 

i A ‘ T nn ViTT* //tt» 7T] 


X 0 T S GF/(F-1) 


(m' -GX 0 T S ) 2 
•exp [ — u s/ 


_ v _ m' -GX 0 T S 

2X 0 T s G 2 F(1+ — — r ±±- ) 


]dm' 


GX 0 T s F/(F-1) 


Since, for m' >GX 0 T g 


-GXqT, m' -GX 0 T Si m , 

X ‘^1\ rn n //m * \ ^ V 1 ' ~ ^ ) 


GX 0 T s F /(F- l) 


GX 0 T S ' GX 0 T S 


then, 


m' -GX ft T c 


2X 0 T,G 2 F(l + - ^:_,r S .' )<2GFm' . 


GXoV/fF-l) 


The exponent in (A4.b) satisfies 
( m' -GX 0 T S ) 2 


m f — GX n T c 

2X 0 T s G 2 F(1+ 2_J_) 

US V I /~1\'T‘T7‘//T7'1W 


> 


GX 0 T s F/(F-1) 


(m' -GX 0 T S ) 2 
2GFm' 


m' 2 -m' GX 0 T S 
2GFm' 


(A4.b) 


m' -GX 0 T S 
2GF 


Also 


then, 


; i + 


m' -GX 0 T, 0/0 m' -GX n T, 

0 s 13/2 v r 0 s 13/2 

X 0 T,GF/(F-1) 1 ^ 1 X 0 T S GF/(F-1) 1 ’ 


e 3 < 


Q-i 


I 


(2^X 0 T s G 2 F) 1 / 2 M' r m ' ~ GX 0 T 3 ,3/2 
1 X 0 T S GF/(F-1) J 


f m' -2X 0 T S G 
‘ ex P [ ^ l dm ' 


(A4.c) 


2GF 
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T * _ 2GF , 

Let u= — - t — - — , so that 

m' -X 0 T 3 G 


2GF 


e 3 


. M' -X 0 T,G _J_ 

< l(Q-l)X 0 T s e XoT>/2F -. /— 1 e" u 
2 V V 7r(F — 1 ) 3 o 


du 


(A4.d) 


r(F-l) 3 o 

The proper value of M' is found by setting the upper limit of the integral in 
(A4.d) such that e 3 is smaller than a desired value. Following are two examples 
with two sets of the parameters used in the experiments described in section III. 

Example 1: 

X 0 T S =1, k eff =0.01, G=500, Q=4. For e 3 <lCr 9 , 

2GF 


M' -X 0 T S G 


< 0.067, orM' >105,000 


Example 2: 

X 0 T S =0.01, k eff =0.01, G=700, Q=4. For e 3 <10" 9 , 
2GF 


M' -X 0 T S G 


< 0.095, or M' > 133,000 
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Figure Captions 


Fig. 1. Block diagram of an optimal PPM receiver. 


Fig. 2. Distribution of number of secondary electrons, m, output by an APD, 
given XT S = the average number of absorbed photons. The solid curves are calcu- 
lated using (7) and the dashed curves correspond to the Gaussian density func- 
tion with mean = GXT S and variance = G 2 FXT S . 


Fig. 3. Probability of bit error vs. average APD gain. The solid curves are calcu- 
lated using (14) and the dashed curves are based on the Gaussian approximation. 
Other parameter values used were: APD load resistor R=1030D, equivalent noise 
temperature T=1100°K, APD surface leakage current Ig=11.9nA and PPM slot 
period T s = 20 ns. 


Fig. 4. Same as Fig. 3, but at a different background radiation level. 


Fig. 5. Probability of bit error vs. average number of photons required to 
transmit one bit of information. The solid curves are calculated using (14) and 
the dashed curves are based on the Gaussian approximation. Other parameters 
were the same as in Fig. 3. 


Fig. 6. Same as Fig. 5, but at a different background radiation level. The Gaus- 
sian approximation used G=450, the optimal value according to Fig. 3. The 
crosses represent the experimentally measured data. 


Fig. 7. Same as Fig. 5, but at a different background radiation level. The Gaus- 
sian approximation used G=450, the optimal value according to Fig. 4. The 
crosses represent the experimentally measured data. 


♦ 













V' 


4> 





Probability of Bit Error 



xT 


n 

r 


4 \ 



